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Abstract: We use a volume-limited sample of quasars in the Sloan Digital Sky Survey (SDSS) DRY 

quasar catalog to identify quasar groups and address their statistical significance. This quasar sample has 
a uniform selection function on the sky and nearly a maximum possible contiguous volume that can be 
drawn from the DR7 catalog. Quasar groups are identified by using the Friend-of-Friend algorithm with 
a set of fixed comoving linking lengths. We find that the richness distribution of the richest 100 quasar 
groups or the size distribution of the largest 100 groups are statistically equivalent with those of randomly- 
distributed points with the same number density and sky coverage when groups are identified with the 
linking length of 70 fe~ 1 M pc. It is sho wn that the large-scale structures like the huge Large Quasar 
Group (U1.27) reported bv I Clowes et al.l (120131) can be found with high probability even if quasars have 
no physical clustering, and does not challenge_the initially homogeneous cosmological models. Our results 
are statistically more reliable than those of lNadathur ( 2013h . where the test was made only for the largest 
quasar group. It is shown that the linking length should be smaller than 50 h -1 Mpc in order for the 
quasar groups identified in the DR7 catalog not to be dominated by associations of quasars grouped by 
chance. We present 20 richest quasar groups identified with the linking length of 70 /i _1 Mpc for further 
analyses. 
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1. Introduction 


The size and mass of large-scale structures 0 (hereafter 
LSS) in the universe have been often used as measures 
of amplitude of density fluctuations on l arge sc ales. 
The CfA Great Wall discovered by Ide Lapparent et ahl 
(Il986h in the CfA 2 survey and the Sloan Great Wall in 


the SDSS surv ey dG ott et a h 2005li c ast doubt on the 


SCDM model dG eller fe Huchralll989t) and the ACDM 
model (ISheth fe Diaferiol 1201 ill , respectively. It has 
been argued that the universe did not have enough time 
for such large structures to form by the present epoch in 
these cosmological models, and that the gravitational 
instability theory for structure formation or the Cosmo¬ 
logical Principle adop ting large-scale homogeneity may 
not b e valid (see also lClowes et al.ll2012tlHorvath et al.l 
1201 3lh 

On the other hand, a strong support for the 
Cosmological Principle comes from the smooth cos- 


CORRESPONDING AUTHOR: H. Song 

1 We use the word ‘large-scale structure’ as a general name 
for structures larger than galaxy clusters. In this terminol¬ 
ogy the large-scale structures of the universe are classified into 
high-density and low-density large-scale structures. The high- 
density large-scale structures are classified into categories of 
superclusters, filaments, and walls depending on their internal 
density and morphology. The low-density large-scale structures 
are voids, and classified into bubbles, and tunnels depending on 
their morphology. In this scheme all the names that have been 
used to describe ‘large’ structures in the universe are unified. 


mic microwave background radiation, which shows 
anisotropies that agrees astonishingly well with those 
of the initiall y hom o geneous isotropic ACDM mode l 
( Clifton et al.l l2012t IPlanck collaboration XVl l2014f) . 
The fluctuation of the galaxy number density on large 
scales is also fo und to agree with the p rediction of the 
ACDM model ( Scrimgeour et, al.ll2012li . 

It is expected that LSS of the universe contain a 
wealth of information on primordial density fluctua¬ 
tions and the history of their growth. This is because 
they represent density fluctuations on linear scales. It 
is then very important to rigorously study the physical 
properties and cosmologic al meaning of observed super¬ 
structures of the universe (lParklll990l : iPark et alJl2012t 
lEinasto et al.l 20TlL 201411. _ 

Recently, 1 Clowes et al l (2013) reported a finding of 
a very l arge quasar group in the SDSS DR7 quasar 
catalog (I Schneider et alJl2010H with characteristic size 
of ~350 h -1 Mpc, longest dimension of ~870 h _1 Mpc, 
membership of 73 quasars, and mean redshift z = 1.27. 
They claimed that the quasar group, named U1.27, is 
the largest ‘structure’ currently known in the universe. 
They argued that their finding might raise a question 
on validity of the large-scale homogeneity assumption, 
which is, in the form of Cosmological Principle, one of 
the fundamental assumptions of the standard cosmol¬ 
ogy. Their claim falls on the same kind of view applied 
to CfA and Sloan Great Walls. 
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However, iNadathui] (120131) pointed out that the sta¬ 
tistical significance of the h uge q uasa r grou p has not 
been properly studied in lClowes et al.l ( 2013 ). and that 
Clowes et al.’s claim on the cosmological implication 
of existence of the large quasar group needs to be ex¬ 
amined quantitatively. Nadathur showed that the al¬ 
gorithm used to identify the quasar groups frequently 
finds such large-size structures even in homogeneous 
simulations of a Poisson point process with the same 
number density as the quasar catalog. They concluded 
that Clowes et al.’s interpretation of U1.27 as a physical 
‘structure’ is misleading. 

Even though it is cosmologically in teresti ng 
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: i Clowes & Gan 

lousanol 1991) Komberg et al. 

1996 

: Williger et al. 

200i Clowes et al. 20ia 2013}; 


terpreting the quasar groups identified with specific 
criteria. This is because the size and richness of LSS 
depend sensitively on how they are identified. Evolu¬ 
tion of the objects used to trace LSS and evolution of 
LSS themselves can also complicate the interpretation 
of identified LSS. In particular, the quasar phe- 


activity can affect the results ( 

Kombere; et al. 199^ 

Wold et al. 20001: Miller et al. 

2004} Sochting et al. 

2004 

Coldwell & Lambas 2006: 

Hutchings et al. 

2009 

: Lietzen et al. 
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, 2011; 

Bradshaw et al. 

2011 

: Portinari et al. 

2012 
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krumne et all 2013: Fanidakis et all 2013; Shen et al. 

2013} DiPonmeo et al. 2014[ Karhunen et al. 2014 

Song et al. 2015h - Quasar samples usually have very 


low spatial number density, which can introduce shot 
noise effects in group identification. This situation 
demands the use of well-defined samples and rigorous 
statistical tests. 

In this paper we identify quasar groups using a 
volume-limited sample drawn from the SDSS DR7 
quasar catalog that has the maximum contiguous an¬ 
gular size. We present the statistical significance of 
not only the largest quasar group, but also the quasar 
groups within a large range of richness. We also address 
a question whether or not the huge quasar group found 
by Clowes et al. is a statistically-significant object in 
terms of multiplicity and size. We find the critical link¬ 
ing length (LL hereafter) for quasar group identification 
which can result in a catalog dominated by physically 
clustered groups. 


2. SDSS Quasar Catalog 


We use the fifth edition of the SDSS quasar catalog from 
SDSS DR7, which is a compilation of q uasars pro bed in 
the SDSS-I and SDSS-II quasar survey (Schneider et al.l 


2 In this paper we call the large loose groupings of quasars 
‘quasar groups.’ This terminology can give a misleading idea 
that they are com pact and physically i nteracting internally as 
in galaxy groups. lEinasto et al.l <120141) used ‘quasar systems’ 
to avoid the confusion. The name ‘quasar association’ may also 
be an appropriate name for these loose quasar groups since the 
name ‘stellar association’ has been already used in astronomy 
for loose groups of stars that are not gravitationally bound. 


l20ld) . The catalog contains 105783 spectroscopically 
confirmed quasars in the sky of area covering approxi¬ 
mately 9380deg 2 . They are brighter than Mi = —22.0 
(corrected for the Galactic extinction, and RT-corrected 
to z = 2 in a cosmology with H 0 = 70km/s/Mpc, 
D m = 0.3, and Da = 0.7) and have at least one broad 
emission line with full width at half maximum larger 
than lOOOkm/s or interesting/complex absorption fea¬ 
tures. The quasars show a wide redsliift distribution 
from 0.065 to 5.46, but the majority of quasars (76%) 
is concentrated below redsliift 2. They have the appar¬ 
ent i magnitude of 14.86 < i < 22.36, where the lower 
limit comes from the maximum brightness limit of the 
target selection on quasar candidates to avoid satura¬ 
tion and cross-talk in the spectra. The catalog does 
not include several classes of active galactic nuclei such 
as Type 2 quasars, Seyfert galaxies, and BL Lacertae 
objects. 

Even though the catalog contains a large number 
of quasars, not all the quasars are suitable for sta¬ 
tistical analyses for several reasons. It is mainly be¬ 
cause the quasar target selection for the follow-up spec- 
troscopy was not done uniformly in its early version (see 
Richardset ah| |2006l : IShen et al'll2007l for more details). 
Shen et al.l (|201lh provided a parameter, in the 10 th col¬ 
umn in their catalog, named as ‘uniform flag’ to take 
the uniformity problem into account. The catalog of 
Schneider et al. also has its own uniform flag in the 
35th column. The uniform flag is assigned based on 
the achievement of the target selection algorithm. But 
Shen et al. made more careful evaluation of the selec¬ 
tion effects by finding out quasars that are uniformly 
selected as high-z candidates but later identified as low- 
z quasars. Shen et al. has one more class than Schnei¬ 
der et al. as follows: 0 = not in the uniform sample; 
1 = uniformly selected and with the galactic extinction- 
corrected i magnitude less than 19.1 at z < 2.9 or 20.2 
at z > 2.9; 2 = uniformly selected by the QSO_HiZ 
branch in the algorithm and with measured spectro¬ 
scopic redsh ift z < 2.9 and i > 19.1. By selecting those 
whose IShen et akl (|201lh ’s uniform flag is 1, we con¬ 
struct a uniform sample of 59679 quasars with Galactic 
extinction corrected i magnitude less than 19.1. It re¬ 
jects about a half of the quasars in the original catalog, 
mostly at z < 2.9. 

In F igure 1 we plot the quasars in the lSchneider et ahl 
(f2010h ’s catalog in black. Only the northern hemisphere 
of the SDSS survey regions is adopted and the SDSS 
survey coordinates rj and A are used for the plot as 
the survey boundaries are more simply defined in this 
coordinate system. The red points in the figure are 
quasars marked as those of the uniform subsample by 
Shen et al. 

Since the sky distribution of the uniform quasar sam¬ 
ple is not perfectly contiguous, we define a new bound¬ 
ary for our sample by excluding isolated patches and 
trimming the jagged boundaries to make one clean con¬ 
tiguous region. The green line in Figure 1 defines the 
region under our study. The area of the region is 1.754 
steradian. This is nearly the maximum possible con- 
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tiguous sample can be drawn from the DR7 catalog. 

Figure 2 shows the uniform-sample quasars located 
within our mask in the i-band absolute magnitude and 
redshift plane. In Figure 3 it can be seen that the co¬ 
moving space number density of quasars is roughly uni¬ 
form between z = 0.8 and 1.8, but is significantly higher 
b etween z = 0. 1 and 0.8. 

ISong et all (1201511 has compared the distribution of 
the SDSS quasars at 0.4 < z < 0.7 with the smooth 
density field of the galaxies in the SDSS DR12 BOSS 
survey and found that the probability of finding a 
quasar increases monotonically as the local galaxy num¬ 
ber density increases with no significant redshift depen¬ 
dence. This implies that LSS found in the distribution 
of quasars should be similar to those in the galaxy dis¬ 
tribution. Then a quasar sample having a constant co- 
rnoving number density should be as useful as a volume- 
limited galaxy sample in finding LSS. 

In order to make a volume-limited quasar sample 
with uniform density we apply an absolute magnitude 
cut line given by 

= 0.7216s 2 - 3.505s - 22.278 (1) 


in the range of 0.1 < s < 0.8. The cut is shown as a blue 
line in Figure 2. The resulting subsample has a nearly 
uniform number density of 1.31 x 10 _6 (/i _1 Mpc) _3 in 
a wide range of 0.1 < s < 1.8 as shown by the red his¬ 
togram in Figure 3. The corresponding mean quasar 
separation is 91.5 h -1 Mpc. By varying the absolute 
magnitude cut one can effectively correct a quasar sam¬ 
ple for the luminosity evolution. The use of uniform 
comoving density sample makes us to sample the same 
kind of cosmological objects over a wide redshift range. 

Our final sample contains 32276 quasars with i < 
19.1 in the redshift range 0.1 < z < 1.8 co vering 
5758 deg 2 . On the o ther hand, IClowes et al.l (l2013lf 
and iNadathurl ( 20131 1 did not use the sample made 
to meet the precondition of uniformity for statistical 
works. They just applie d an i-ban d apparent magni¬ 
tude cut of i < 19.1 to lSchneider et ah (120101') ’s sample, 
where i here is not corrected for the Galactic extinction. 
This made their sample affected b y non-uniform quasar 
target selection which IShen et all (1201 i ll corrected, and 
the structures identified are affected by the Galactic 
extinction. Our sample is from the homogeneous sam¬ 
ple and has a true flux limit. It is also much larger 
in the sky area (5758 deg 2 versus 3725 cleg 2 ) and in 
redshift range (0.1 < z < 1.8 versus 1.0 < z < 1.8). It 
makes our s ample contain 1.72 t i mes more quasars than 
that used bvIClowes et al.1 (|2013ll . INadathurl (12013H . and 
lEinasto et ah ( 201411 . 


3. Identification of Quasar Groups 

To identify groups of quasars we apply the Friend- 
of-Friend (FoF) algorithm to the ‘uniform’ constant- 
number density sample described in the previous sec¬ 
tion. We hrst use LL= 70h _1 Mpc for a comparison 
with IClowes et al.l (l2013H ’s results. This corresponds 
to the threshold overdensity of <5 = ( d/LL ) 3 — 1 = 



Fig ure 1. Distribu t ion o f all the quasars in the origi¬ 
nal [Schneider et ah[ j[2010|Vs cata log (black points) and the 
quasars in the [Shenetahl (120111 1 ’s uniform sample (red 
points) in the SDSS survey coordinates. A contiguous re¬ 
gion defined by the green line is the region under analysis 
in our study. 



Figure 2. Distribution of the uniform-sample quasars in the 
redshift versus i-band absolute magnitude plane. The blue 
line is a low-redshift absolute magnitude cut used in our 
study making the sample nearly constant in comoving num¬ 
ber density between z = 0.1 and 1.8. 
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Figure 3. Comoving space number density of quasars. The 
black line is for the whole uniform sample, and the red line 
is when the low-redshift absolute magnitude cut (Eq. 1) is 
applied. The blue dashed line is the mean density between 
z = 0.1 and 1.8. 


(91.5/70) 3 — 1 = 1.23, where d, is the mean quasar 
separation. We find 3233 quasar groups having 3 or 
more members. Figure 4 shows the positions of top ten 
richest groups on the sky (upper panel) and in the x-y 
plane of Equatorial coordinate frame (bottom panel). 
The second richest one labelled with ‘2’ has 70 member 
quasars and the maximum extent of 703 h _1 Mpc. The 
maximum extent is defined as the comoving distance 
between the most distant two members in a group. This 
group corresponds to huge Large Quasar Group, U1.27, 
reported by Clowes et al. The group identified by us has 
3 new members and does not have 6 members compared 
to Clowes et al.’s list. We discover a rich quasar group, 
marked with G’ in Fig. 4, that contains 95 members 
and the maximum extent of 712 h~ 1 Mpc. This group 
was not found by previous studies because its center is 
located at z = 0.366, which is outside the redshift limit 
of 1.0 < x <1.8 adopted by Cl owes et al., I Nad at hurl 
(l2013f) . and lEinasto et al.l (l2014fi . This group is much 
richer and has a size larger than U1.27. A catalog of 
twenty richest quasar groups is given in Table 1 for fu¬ 
ture analyses. 

To check if our group finding is consistent with previ¬ 
ous works we perf orm the FoF group find ing to a sample 
from the original [Schneider et al.l (|2010l Vs catalog that 
is limited only by i < 19.1 not corrected for Galactic 
extinction and a redshift range of 1.0 < z < 1.8. In this 
case the richest group is found to have 73 members that 
are exactly the same as the member quasars of U1.27 
listed in Table 1 of Clowes et al. This proves that the 
difference in our quasar catalog is not due to the group 
finding algorithm but to the sample definition. We also 
compared our quasar groups with those identified by 


Table 1 

The twenty richest SDSS DR7 quasar groups identified 
with the linking length of 70 /i _1 Mpc or overdensity 
threshold of 1.23. 


ID 

a(deg) a 

5(deg) a 

z a 

r a 

N b 

L c 

1 

170.08 

29.33 

0.366 

1012 

95 

712 

2 

163.53 

14.47 

1.244 

2757 

70 

703 

3 

196.45 

39.95 

1.112 

2546 

64 

536 

4 

195.98 

27.09 

1.554 

3198 

62 

625 

5 

166.91 

33.88 

1.075 

2485 

60 

570 

6 

156.69 

38.63 

0.580 

1519 

59 

536 

7 

219.41 

27.68 

0.598 

1560 

56 

459 

8 

204.40 

13.04 

1.215 

2712 

55 

582 

9 

187.60 

43.55 

1.404 

2995 

53 

651 

10 

232.82 

23.19 

1.511 

3141 

53 

556 

11 

219.69 

19.68 

0.764 

1909 

53 

640 

12 

226.57 

16.83 

1.062 

2463 

51 

470 

13 

141.58 

47.90 

1.157 

2620 

50 

526 

14 

208.63 

25.81 

1.236 

2745 

48 

446 

15 

137.82 

21.70 

0.957 

2279 

47 

414 

16 

212.06 

28.23 

1.080 

2493 

47 

754 

17 

127.45 

20.13 

1.351 

2917 

46 

545 

18 

157.40 

22.01 

1.493 

3117 

45 

576 

19 

158.17 

15.82 

0.744 

1869 

45 

456 

20 

231.00 

47.84 

1.508 

3138 

44 

459 


The full catalog with information on members of groups can be 
downloaded at http://astro.kias.re.kr/Quasar/DR7_20groups/ 
a Coordinates of the geometric center of groups. The mean co¬ 
moving distance f is in units of h~ ■’^Mpc. 
b The number of member quasars. 
c The maximum extent of groups in h~ 1 Mpc. 


lEinasto et al.1 ( 2014 see their Table A.3) who also used 
Schneider et al.’s catalog constrained by i < 19.1, and 
found a very good agreement at z > 1.0. We choose our 
sample instead of that used by previous works because 
the observed sample must be uniformly defined for sta¬ 
tistical hypothesis tests against random catalogs. The 
observed quasar distribution has random perturbations 
along the line-of-sight due to the random redshift er¬ 
ror which is t ypically about 0.004 in the case of SDSS 
DR7 catalog (ISchneider et alJf2010ll . Therefore, it is 
necessary to examine if the quasar group finding is sig¬ 
nificantly affected by this error. We will discuss it in 
the next section. 


4. Comparison with random distributions 

We generate random point sets which have the same 
number density and the same angular coverage on 
the sky as our quasar sample. The random catalogs 
also have the mean point separation of 91.5h -1 Mpc 
throughout the entire redshift range from z = 0.1 to 1.8. 
We make 10000 random mock samples. For each mock 
sample we find groups of points in the same way as for 
the quasar sample. We then measure richness and size 
of the groups in quasar sample and random mock sam¬ 
ples. Richness is the number of member quasars/points 
belonging to a given group. Size is the maximum co¬ 
moving extent of the group. 

In the upper panel of Figure 5 we plot the cumulative 
distribution function (CDF) of richness. It is the num- 
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Figure 4. (upper panel) Ten richest quasar groups detected 
with the linking length of 70hr 1 Mpc. The second rich¬ 
est one in magenta on the bottom right labelled with 2 
correspon ds to the huge Large Quasar Group, U1.27, of 
IClowes et al.l (|2013f l. (bottom panel) Same as the upper 
plot but in the Cartesian coordinates defined in the Equa¬ 
torial coordinate system. Grey points are all the quasars in 
our sample. 



N . 



Figure 5. Cumulative distribution functions of richness (up¬ 
per panel) and maximum extent (bottom panel) of the 
groups identified in the observed SDSS quasar catalog (red 
circles) and in random point catalogs (solid line with error 
bars). A set of 10000 random point catalogs is used for the 
comparison. The group correspondin g to the liuge_ Large 
Quasar Group, U1.27, reported bv IClowes et all (2013) is 
marked with downward arrows. 


ber of groups whose richness is greater than or equal to 
N g . The circles are from the observed quasar sample. 
Richness of the group corresponding to U1.27 group is 
70, and is marked with an arrow. The solid line and 
error bars are the mean and standard deviations ob¬ 
tained from 10000 random catalogs. Similarly, the bot¬ 
tom panel shows the number of observed quasar groups 
whose size is greater than or equal to L, and the results 
from 10000 mock catalogs. We find that U1.27 is 703 
/i~ 1 Mpc long (marked with an arrow) and is only the 
fourth largest one. It can be seen that the richness and 
size distributions of the observed quasars are statisti¬ 
cally indistinguishable with those of random catalogs. 

To estimate the probability of finding groups like the 
observed quasar groups in random catalogs we perform 
a x 2 test. We define 

X 2 = ~ <k) 2 /°f ( 2 ) 

»=l 

as a measure of the deviation of a particular CDF <f> 
from the average CDF cj) of random catalogs. Here cq 
is the standard deviation of the CDF’s of random cat¬ 
alogs in the i-th bin, and the summation is over the 
range where the average CDF of random catalogs is 
between 3 and 100. The choice ensures that a large 
part of CDF is compared. The histograms in Figure 
6 are the distributions of x 2 of 10000 random catalogs 
and the arrows correspond to the observed SDSS quasar 
catalog generated with LL= 70/i _1 Mpc. The probabil¬ 
ity of having x 2 greater than x 2 bs i n random catalogs 
is 58.1% and 92.4% for richness and size distributions, 
respectively. Therefore, it can be concluded that the 
quasar groups identified with LL= 70/i _1 Mpc are con¬ 
sistent with those of random point distributions. 

On the other hand, we find that the probability of 
finding a random group whose richness is equal to or 
greater than 70 is 92.3%. Similarly, the probability of 
finding a random group with the maximum extent equal 
to or greater than 703 /i -1 Mpc is 100%. Therefore, 
both richness and size of the objects like U1.27 are of 
no surprise even in random distributions of points with 
no physical clustering when LL is as large as 70ft. -1 Mpc. 

Since quasars do have intrinsic spatial clustering, it 
is expected that these distribution functions are sta¬ 
tistically different from those of random catalogs for 
some smaller LL or for higher density thresholds. Fig¬ 
ure 7 shows the CDF’s of richness and size of quasar 
groups when LL is taken to be 60, 50, and 40 /i _1 Mpc. 
Also shown are the average CDF’s (solid lines with er¬ 
ror bars) of random point groups identified with the 
corresponding LL. It is clear that the CDF’s of ob¬ 
served group richness and size become more different 
from those of random catalogs as LL decreases. 

It will be then interesting to search for the critical 
LL for which the observed CDF becomes significantly 
different from those of random points. We calculate the 
probability that the CDF of random point groups hap¬ 
pens to show a deviation from the mean CDF as large 
as or greater than the difference of the observed CDF 
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Figure 6. Distributions of \ 2 defined as in Eq. (2) for 
random point samples. The upper panel is for richness of 
groups identified with the linking length of 70 /) _1 Mpc, and 
the bottom panel is for the maximum extent of groups. The 
\ 2 for the observed sample are marked with arrows. 

from the mean. This is estimated by counting the cases 
with Xran — X'obs' A high probability means that the 
observed groups are statistically consistent with ran¬ 
dom ones. Figure 8 shows the probability as a func¬ 
tion of the LL used to find groups. For both richness 
and size it can be seen that the observed groups are 
consistent with those of random catalogs when LL is 
chosen to be larger than about 50h _1 Mpc. Clowes et 
al. chose LL= 70/i _1 Mpc and it is now clear that the 
choice makes the group catalog dominated by randomly 
connected groups. Figures 7 and 8 indicate that the 
groups identified with LL< 50h -1 Mpc are likely to be 
physically-clustered genuine groups of quasars. 

As mentioned in the previous section the random 
redshift error of order of 0.004 seems problematic in 
finding physical structures as the error corresponds 
to ~ 12h~ 1 Mp c . We use the redshifts calculated by 
iHewett fe Wildl (|2010l hereafter HW10) whose redshift 
errors are much smaller and redo the group findings. 
According to the error estimates for the redshifts of 
HW10, the RMS redshift error for quasars at 0.1 < z < 
1.8 is found to be only 0.00024, more than 10 times 
smaller than the nominal value for the SDSS DR7. 

It will be a good consist ency check to compare two 
group catalogs derived from lSchneider et~ahl (|2010d and 
HWlO’s redshifts. We find that the new quasar group 
cat alog is very close to t he previous one obtained us¬ 
ing |Schneider]eFjjI] (l2010h ’s redshifts. All major groups 
are detected again and only a few members are affected. 
As a result, the richness and size distributions change 
only slightly. It is reassuring that the quasar groups 
and the statistics are robust against the redshift error 


Figure 7. Cumulative distribution functions of richness (up¬ 
per panel) and size (bottom panel) when the linking length 
for group identification is taken to be 60, 50, and 40 
/i -1 Mpcs. Circles are from observation, and solid lines with 
error bars are from 1000 random catalogs. 

of lSchneider et al.1 (|2010t h When LL = 50/i _1 Mpc, we 
still find a very good agreement of both distribution 
functions. Therefore, all of our results are essentially 
unchan ged by using mor e accur ate quasar redshifts. We 
stick to iSchneider et all d2010l) ’s redshifts in our paper 
since there is practically no change in the conclusions 
and it is easier to compare our groups directly to those 
of previous works. The groups identified in the red¬ 
shift space are also affected by the peculiar velocities of 
quasars. Influence of peculiar velocity on group iden¬ 
tification is to some extent unavoidable. But this is 
relatively small because the RMS pairwise peculiar ve¬ 
locity of galaxies is only about 300 km/s. 

5. Summary and Conclusions 

We make a wide contiguous volume-limited sample of 
quasars statistically uniformly selected from the SDSS 
DR7 quasar catalog and identify quasar groups using 
the FoF algorithm. The statistical properties of the ob¬ 
served quasar groups are compared with those of ran¬ 
dom point sets. We make the following findings. 

1. The observed SDSS DR7 quasar groups identi¬ 
fied with LL= 70/i _1 Mpc are statistically consistent, in 
terms of size and richness, with those of random points 
clustered by chance. The size or richness distribution 
functions of those quasar groups can be found in ran¬ 
dom catalogs with 92.4% and 58.1% probabilities, re¬ 
spectively. 

2. The probabilities of finding groups larger or 
richer than the large quasar group U1.27, reported by 
IClowes et al.l (120131) . in a random point catalog are 
100% and 92%, respectively. 
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Figure 8. Probability of finding random catalogs whose dis¬ 
tribution function deviates from their mean more than or 
equal to that of the observed sample. The deviation is mea¬ 
sured by x 2 defined in Eq. (2). Probability is measured 
for richness (solid line) and size (dashed line) distribution 
functions as a function the linking length. 


There have also been several studies trying to 
define an d measure t h e homogeneity scale o f the 
universe ( Yadavet_aLl l20ld : iMarinoni et, al.l l2012t 
IScrimgeour et alJl2012|) . However, definition of the ho¬ 
mogeneity scale has been either arbitary or dependent 
on observational samples. We point out that, in the 
current popular cosmological models that adopt the 
Cosmological Principle, the notion of the critical scale 
of homogeneity does not make sense. Homogeneity is 
only achieved asymptotically in these models as the 
scale of observation increases, and there is no charac¬ 
teristic scale intrinsic to the universe above which the 
univ erse can be suddenly regarded homogeneous (see 
also lNadathuill2013l ). Instead of looking for arbitrarily- 
defined homogeneity scales or naively examining the 
size of the largest structures the observed LSS should 
be statistically compared with the structures identi¬ 
fied in the same way in the mock survey samples ob¬ 
tained from simulated initially-homogeneous universes. 
In our forthcoming paper we will use mock quasars 
formed in simulated universes to identify large groups 
of quasars and statistically compared their properties 
with those of observed quasars. We will examine if the 
large structures can be explained in the framework of 
the standard hierarchical structure formation model in 
the given Hubble time. 


3. The quasar groups are statistically dominated by 
physically clustered quasars when LL for group identi¬ 
fication is less than 50 /i” 1 Mpc. 

We emphasize that the observed LSS should not be 
used to draw cosmological conclusions without making 
elaborate statistical tests. Size and richness of cosmic 
structures depend sensitively on how those structures 
are identified. For example, one can always find struc¬ 
tures as large as the survey size if the threshold density 
level used to identify them is lowered close to the per¬ 
colation density. Even in a random point distribution 
without any physical clustering an infinitely long struc¬ 
ture can be identified when LL is increased to 0.864 
times the mean particle separation, which is the critical 
value for percolation. It is conceivable that one would 
claim discovery of groups of galaxies or quasars even 
larger than the Sloan Great Wall in existing or future 
redsliift surveys. The lesson of the past is that one will 
again need to make statistical tests before drawing con¬ 
clusions on cosmological implication of those objects. 

Our results suggest that the large quasar groups 
identified with LL > 50/i -1 Mpc can be considered 
as associations of smaller physically-clustered quasar 
systems that are connected by a very large linking 
length. The smaller quasar systems are similar to the 
supercluster complexes in the lo c al universe, as was 
sugges ted in lEinasto et al.l (12014P . iHutsemekers et all 
(2014) showed that quasar spin axes in U1.27 seem to 
be aligned parallel to this structure. However, from 
their Figure 6 one can see that alignment is better if 
this quasar group is divided into smaller systems, in 
agreement with our interpretation of quasar systems. 
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